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The JAK-STAT pathway represents a finely tuned orchestra 
capable of rapidly facilitating an exquisite symphony of 
responses from a complex array of extracellular signals. This 
review explores the evolution of the JAK-STAT pathway: the 
origins of the individual domains from which it is constructed, 
the formation of individual components from these basic 
building blocks, the assembly of the components into a 
functional pathway, and the subsequent reiteration of this 
basic template to fulfill a variety of roles downstream of 
cytokine receptors. 



Introduction 

The ability to respond to external signals is essential for life. At 
the level of the whole organism these signals are principally cap- 
tured by peripheral sensory organs, which transmit them via the 
nervous system to the brain to mediate an appropriate response. 
At the cellular level, these signals are generally received by a myr- 
iad of cell-surface receptors that utilize a range of intracellular 
signaling pathways to communicate with the nucleus. This rap- 
idly impacts on target gene transcription and mediates the desired 
cellular outcome, such as proliferation, differentiation, activation, 
and survival. The Janus kinase-signal transducer and activator 
of transcription (JAK-STAT) pathway represents a major mecha- 
nism used to transmit signals from extracellular receptors to the 
nucleus. The canonical pathway consists of JAK and STAT pro- 
teins activated by receptor ligation and inactivated by negative 
regulators, including SH2-containing protein tyrosine phospha- 
tase (SHP) and suppressor of cytokine signaling (SOCS) proteins. 

Each component of the JAK-STAT pathway consists of a col- 
lection of protein domains that collectively provide their distinct 
functional properties. Moreover, the pathway and its components 
have been replicated many times over, comprising four JAKs, 
seven STATs, eight SOCSs and two SHPs downstream of over 50 
cytokine receptors in mammals.' '' This review aims to capture in 
broad brush strokes the current state of knowledge on how such 
a complex signaling system has developed. In particular it fol- 
lows the evolution of each individual component, discusses their 
assembly into a functional system and examines the subsequent 
multiple reiterations of the basic pathway. 
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JAK-STAT Pathway Components 

To understand the evolution of the various components of the 
canonical JAK-STAT pathway, it is necessary to be familiar with 
their constituent protein domains, and how these modules are 
assembled into individual components (Fig. 1). 

At the core of the JAK-STAT pathway is a binary code, based 
on specific tyrosine residues that are either unphosphorylated (off) 
or phosphorylated (on). Three of the protein domains utilized 
in the JAK-STAT pathway are specifically involved in the writ- 
ing, reading and erasing of this code: the protein tyrosine kinase 
(PTK), Src-homology 2 (SH2) and protein tyrosine phosphatase 
(PTP) domains, respectively (Fig. lA). The PTK domain con- 
sists of -250-300 residues, including an ATP-binding site and 
a catalytic region, and is able to phosphorylate tyrosine residues 
with a specificity determined by access to substrate and sequence 
context, thereby writing the code.^ The SH2 domain is a -100 
residue motif that binds to phosphotyrosine residues in a man- 
ner dependent on the specific residues around the tyrosine espe- 
cially those downstream,* providing a mechanism for reading the 
code.' Finally, the PTP domain consists of -280 residues, with a 
signature CX^R motif, and is capable of dephosphorylating phos- 
photyrosine residues with which it is able to effectively interact, 
thereby erasing the code.'"'" 

JAK-STAT signaling components also consist of several other 
domains that provide accessory functions. The PERM domain, 
named after the founding proteins identified with this domain 
(band 4.1, ezrin, radixin and moesin), consists of three sub- 
domains, F1-F3, that are structurally similar to ubiquitin, acyl- 
CoA binding and plecstrin homology-phosphotyrosine binding 
domains, respectively.'^ The FERM domain is able to mediate 
protein-protein interactions, including adaptor and scaffolding 
interactions with membrane bound proteins." FERM domain 
containing proteins can be divided into three broad groups: (1) tal- 
ins and kinesins, (2) ERMS, GEE, kinases and phosphatases and 
(3) myosins and KIRT.''' The SOCS box is a -40 residue motif 
that mediates interactions with proteasomal degradation pathway 
components — particularly ubiquitin ligases — thereby regulating 
protein half-life.' The DNA-binding domain present in STAT pro- 
teins possesses an S-type immunoglobulin fold seen in members 
of the "p53 clan" — that includes p53 and RUNT — which medi- 
ates sequence-specific binding,'''''^ while the coiled-coil domain is 
involved in protein-protein interactions. Finally, the transactivation 
domain (TAD) — of which there are several varieties — is impor- 
tant for co-factor recruitment and transcriptional responses.'^ 

These domains are assembled to form the individual compo- 
nents of the JAK-STAT pathway (Fig. IB): 
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Figure 1. Architecture of JAK-STAT pathway components. (A) The constituent domains of JAK-STAT pathway components and their functionality. The 
PTK, SH2 and PTP domains contribute to the write, read and erase functionalities of the JAK-STAT pathway, as indicated, while the others have acces- 
sory roles. (B) The JAK-STAT pathway components and their domain assemblage. 



(1) JAKs: Possess N-terminal FERM and SH2 domains, 
which mediate interactions with cytokine receptors and other 
signahng proteins, respectively.'^" A centrally located "pseudo- 
kinase" domain — so-called due to its homology to PTK domains 
but (now disputed"") absence of catalytic function — plays a regu- 
latory role,'^ while a classical PTK domain is positioned at the 
C-terminus.'^'^' 

(2) STATs: Consist of a unique N-terminus, with roles in 
nuclear translocation and protein interactions,"* followed by 
coiled-coil, DNA-binding, SH2 and variable TAD domains.'* 
The coiled-coil domain is involved with nuclear export and regu- 
lation of tyrosine phosphorylation, the DNA-binding domain 
mediates recognition of sequences related to TTCN3_^GAA in 
the promoters of responsive genes, the SH2 domain allows for 
specific recognition and docking to phosphotyrosines on cyto- 
kine receptors, JAKs and other STAT molecules,'" while the 
divergent C-terminal TAD mediates transactivation.^" 

(3) SHPs: Possess tandem N-terminal SH2 domains that 
bind specifically to key substrates, including cytokine receptors, 



a catalytic PTP domain and a divergent C-terminal region that 
contains several tyrosine residues that serve as docking sites for 
other signaling proteins when phosphorylated." "' 

(4) SOCSs: Consist of a divergent N-terminal region of 
largely unknown function, a central SH2 domain responsible 
for binding to cytokine receptors, JAKs and other substrates 
and a C-terminal SOCS box domain that facilitates protein 
degradation of the cytokine receptor signaling molecules.^* 
Specific SOCS members also possess additional motifs, such 
as the kinase inhibitory region (KIR) of SOCSl and SOCS3, 
and the N-terminal conserved region (NTCR) of SOCS4 and 
SOCS5.'"'' 

The JAK-STAT Pathway 

The individual components interact functionally in the canoni- 
cal JAK-STAT pathway (Fig. 2). This pathway is activated by 
cytokine receptors in response to extracellular ligand binding, 
which leads to conformational changes in the receptor that 
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Figure 2. The canonical JAK-STAT pathway. (A) Ligand binding to a cytol<ine receptor results in dimerization and conformational changes that causes 
activation of JAKs associated with the Box 1 and Box 2 domains of the receptor intracellular domain. (B) The activated JAK is able to phosphorylate 
tyrosines on the receptor complex to create docking sites for signaling proteins, such as STATs. These are in turn phosphorylated, and then dimer- 
ize and translocate to the nucleus, where they bind to specific sequences in the promoters of responsive genes. (C) The pathway is extinguished by 
pre-existing SHPs, which serve to dephosphorylate tyrosine residues, or induced SOCSs that can act via several mechanisms, including JAK inhibition, 
steric interference at STAT docking sites, or mediating degradation of receptor signaling components. 



activate the intracellularly associated JAKs. These phosphorylate 
the receptor complex to recruit STAT proteins, which are in turn 
phosphorylated, allowing dimerization and subsequent translo- 
cation to the nucleus to influence the transcription of relevant 
genes. The activated JAK-STAT pathway is then extinguished 
by negative regulators recruited to the receptor signaling com- 
plex, including latent SHPs or induced SOCSs. Additionally, the 
individual components are also able to signal via alternate "non- 
canonical" mechanisms. ^^'^^ 

Building the Building Blocks 

The first obligate step in the evolution of the JAK-STAT path- 
way was the generation of each of the individual domains that 
make up the various components of the pathway. Only 136 pro- 
tein domains are metazoan-specific, with the majority of protein 
domains (> 400) emerging early in eukaryote evolution,^' which 
is also the case for the constituent domains of the JAK-STAT 
pathway (Fig. 3). 

PTK domain. There is genomic evidence that PTK domains 
emerged in early eukaryotes or indeed prior to eukaryotes, as 
a modification of pre-existing kinase domains with alternate 
amino acid specificity.'^'^^" '^ PTK domains are certainly present 
in a broad range of extant non-metazoans, such as amoebozoa 



(such as the slime mold, Dictyostelium discoideum) , basal holozoa 
(such as the filaseterea, Capsaspora owczarzaki) and choanofla- 
gellatesa (e.g., Monosiga brevicollis), and have become extremely 
widespread in the metazoan lineage due to the myriad of path- 
ways that exploit tyrosine phosphorylation.^^ 

SH2 domain. This domain, present in all JAK-STAT path- 
way components, can be traced to the origins of eukaryotes. 
Proto-SH2 domains are found in present-day plants (e.g., 
Arabidopsis thaliana) and slime molds, with the number of SH2 
domains containing proteins expanding to > 100 in higher verte- 
brates.'*'"''^'' Interestingly, the SH2 domains of the JAK and STAT 
families are closely related when compared with those of other 
SH2 domain-containing proteins.^'^^ 

PTP domain. The PTP domain is present in prokaryotes, bac- 
teria and archaea, and also basal eukaryotes such as malaria (e.g., 
Plasmodium berghei), and indeed may have also been captured 
by several viruses."^ Like the PTK domain, this domain seems 
to have evolved by modification of a pre-existing enzyme to one 
with altered specificity for phosphotyrosines. Curiously, there are 
relatively few examples of PTP domains in extant representa- 
tives of plants and slime molds."'"* However, the PTP domains 
have expanded considerably since the divergence of holozoa,^' 
again reflecting the increasing importance of phosphotyrosine- 
mediated signaling along this lineage. 
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Figure 3. Assembly of tiie JAK-STAT pathiway components during evolution. An abbreviated tree of life for the evolution of eukaryotes (solid black 
lines) including major divergence points (rectangles), timing of WGD events (1R, 2R and 3R) and key organisms (rounded rectangles)/species (shaded 
rectangles) covered by this review. The likely time of derivation of the constituent domains is shown to the left of the tree. The proposed evolutionary 
steps that formed the various JAK-STAT component architectures are shown to the right of the tree, including the subsequent reiteration of individual 
components, the number of which is shown in green. 



FERM. A number of FERM domain containing proteins are 
found in present-day plants and basal eukaryotes,''''"'''"' indicating 
a long evolutionary history. Indeed, there are six FERM domain 
containing proteins in present-day slime molds, although the 
ERMs, GEF, kinases and phosphatase group is absent,''' with a 
total of 22 FERM domain containing proteins in extant insects, 
such as the fruit fly [Drosophila melanogaster) P 

SOCS box. Current genomic evidence suggests that the SOCS 
box is holozoan specific, with a SOCS box containing protein 
seen in present-day choanoflagellates (e.g., Salpingoeca rosetta)}'' 
The number of SOCS box containing proteins has subsequently 
expanded considerably within metazoa.' '' 

Other STAT-specific domains. The individual domains of 
STAT have arisen at different times during evolution. The p53 
clan-type DNA binding domain has been present since early in 
eukaryote evolution,"" with such domains described in extant 
plants,'" early unikonta, such as amoebozoa, and holozoa, includ- 
ing choanoflagellates and ichthyosporea."" Although the roots are 
not as ancient, the coiled-coil domain arose prior to the diver- 
gence of amoebozoa, being evident in present-day slime molds 



and early holozoans. Finally, proto-TAD domains probably arose 
in the bilateria."' 

Component Assembly 

The next important step in JAK-STAT pathway evolution was 
the assembly of these specific domains to generate each indi- 
vidual component. The assembly of domains into new protein 
architectures involves genomic processes such as exon-shuffling, 
exonization, intronization and duplication,*^'''* while sequence 
divergence can modify these domains."' '''' These processes, rather 
than de novo generation of novel domains represents the major 
mechanism for the domain variation of proteins in eukaryotes, 
which is affirmed in the analysis of individual JAK-STAT path- 
way components. 

JAKs. Each of the constituent domains of the JAKs, with 
the exception of the pseudo-kinase domain, was present 
throughout eukaryote evolution. So, how did the unique JAK 
architecture arise? The FERM domain of JAK is most closely 
related to that of focal adhesion kinase (FAK), which consists 
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of an N-terminal FERM domain, SH2-binding domain (not 
actually a SH2 domain), PTK domain and a C-terminal focal 
adhesion targeting (FAT) domain,'* which is evident in basal 
holozoans, such as present-day filasterea (e.g., C. owczarzaki) ?° 
Interestingly, the JAK SH2 domain shows highest homology 
with the N-terminal SH2 domains of two other kinases, SYK 
and ZAP70, which possess dual SH2 domains and an adjacent 
PTK domain.''^ The earliest SYK-like protein is exemplified by 
that in present-day sponge (e.g., Amphimedon queenslandica) ?° 
It remains unclear whether the JAK architecture was generated 
by the addition of an SH2 domain to a FAK-like protein that 
subsequently lost the FAT domain, or the addition of the FERM 
domain to a SYK-like protein, or some fusion of a FAK and 
SYK protein, prior to the divergence of eumetazoa and porif- 
era. Curiously, M. brevicollis contains a protein, UTK05mod, 
which consists of FERM, SH2 and dual kinase domains, that is 
reminiscent of the metazoan JAK proteins, although the kinase 
domain forms a clade with the SYK kinase domain and not 
that the JAK kinase domain.^" However, the similarity of the 
architecture between UTK05mod and JAK may be due to coin- 
cidence as there exist a wide variety of kinase domain contain- 
ing architectures that are specific to M. brevicollis?^ The first 
bona fide JAK-like protein is represented in the basal metazoan 
porifera lineage, typified by the sequence BAA81718.3 found in 
the present-day sponge, Ephydatia fluviatilis, which consists of 
FERM and SH2 domains as well as a single PTK domain related 
to the JAK PTK, although lacking the pseudo-PTK domain.*'' 
The canonical JAK domain architecture inclusive of a pseudo- 
kinase domain appears to be unique within bilateria, with this 
unique topology typified by present-day fruit fly Hopscotch.* 
It is likely that that the pseudo-kinase domain a result of local 
duplication and divergence of the PTK domain within a JAK- 
like sequence, although there exists the possibility that it was by 
modification of the first PTK domain of a UTK05mod-related 
sequence. 

STATs. The STAT family has more ancient roots, with a range 
of STAT-like proteins present in extant eukaryotes. The most 
divergent are the GRAS proteins of the plant kingdom (named 
after the founding members, GAI, RGA and SCR), which possess 
an SH2-like domain that is remarkably similar to the metazoan 
STAT SH2 domain and also a putative DNA binding domain.*' 
By the time of divergence of the amoebozoa, a family of proteins 
can be identified that resemble metazoan STATs, with coiled- 
coil, DNA-binding and SH2 domains, although lacking trans- 
activation domains and the N-terminal sequence, as exemplified 
by those in present-day slime molds. Extant members of the 
cnidara, the sea anemone {Nematostella vectensis) and nematode 
{Caenorhabditis elegans), both contain at least one STAT with 
architectures similar to that seen in slime mold. Furthermore, 
the nematode STAT paralogs also contain a TAD,*''" while the 
sea anemone STAT contains an N-terminal sequence, but not a 
TAD, although it possesses a long C-terminus containing mul- 
tiple tyrosine and serine residues. The canonical STAT, replete 
with both TAD domain and N-terminal sequence, first appeared 
in the bilateria, as exemplified by present-day fruit fly Marelle 
(STAT92E).3'" 



SHPs. The constituent domains of the SHPs were also pres- 
ent throughout eukaryotes evolution, but there is no evidence of 
these domains coming together in the bikonta or amoebozoa. 
However, a hypothetical SHP-like protein, is found in present- 
day choanoflagellates (e.g., XP_001747983 in M. brevicollis), 
consisting of a single SH2 domain followed by a PTP domain,^^'^' 
which suggests that this architecture had evolved by the time of 
holozoa. The canonical SHP was finalised early in metazoan evo- 
lution, with the extant nematode ortholog, ptp-2, and fruit fly 
Corkscrew displaying the unique dual N-terminal SH2 domain 
architecture.'^ This most likely arose via local duplication of the 
adjacent SH2 domain in the precursor as the N-terminal and 
C-terminal SH2 domains of SHP are more similar to each other 
than the SH2 domains of other proteins.^' 

SOCSs. Evolution of the SOCS proteins has followed a 
slightly more unique path. The most ancient SOCS box-contain- 
ing protein architecture consisted of three N-terminal ankyrin 
repeats and a C-terminal SOCS box, similar to that of metazoan 
ankyrin repeat and SOCS box (ASB) proteins." The SOCS SH2 
domain is most closely related to that of dual adaptor of phospho- 
tyrosine and 3-phosphoinositides (DAPPl), which together are 
part of a much larger clade that includes the SHP SH2 domain. 
Presumably domain accretion to fuse a SH2 domain with a 
SOCS box generated the archetypical SOCS. In contrast to other 
components, this prototype underwent several gene duplication 
events that established three SOCS members by the divergence 
of the cnidara. This is exemplified in extant cnidara (e.g., N. vec- 
tensis), which possess three proteins with the signature SOCS 
protein architecture."" These are related to each of the three dis- 
tinct SOCS clades seen in higher vertebrates: SOCS1/2/3/CISH, 
SOCS4/5 and SOCS6/7.* 

When JAK Met STAT and SOCS and SHP: 
The Canonical JAK-STAT Pathway 

At the time of the emergence of bilateria, archetypal members of 
each of the JAK-STAT pathway components had been generated, 
including a single JAK, STAT and SHP, and three SOCS proteins. 
These came together to form the canonical JAK-STAT signal- 
ing pathway prior to the divergence of protostomia, as exempli- 
fied by the present day fruit fly Hopscotch-Marelle- Corkscrew/ 
Socs36E pathway downstream of Domeless receptors. Several 
of the key innovations that were required to bring the various 
components together were, first, the association of the prototypal 
JAK with a membrane-bound receptor that lacked intracellular 
tyrosine kinase activity, second, appropriate tyrosines with which 
to recruit STATs, SHPs and SOCSs, and third, the incorporation 
of a negative feedback loop via the SOCS family. These innova- 
tions provided the ability to both recruit and activate the writer, 
the "paper" on which to write and a vehicle to inducibly block the 
writer or reader and/or act as a "paper shredder" to extinguish 
the signal, thereby complementing the erasers. Surprisingly little 
was required to establish these innovations. A few mutations to 
generate a cluster of intracellular prolines would have been suffi- 
cient for a cytokine receptor to be able to interact with the FERM 
domain of a prototypal JAK.'^''^ A similar process could generate 
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tyrosines in an appropriate sequence context, and a small num- 
ber of mutations in the promoter region of a pre-existing SOCS 
gene would have been enough to generate an appropriate STAT 
binding site (TTCNj_^GAA) to enable their recruitment into 
an inducible feedback loop. Of course, each of the components 
would have necessarily participated in other processes prior to 
their recruitment into the canonical JAK-STAT pathway. Indeed, 
the legacy of this is evident in some of the non-canonical roles 
played by the pathway components in alternate signaling para- 
digms, such as downstream of receptor tyrosine kinases. ^''■^'''^•'''''^ 

Making the Most of a Good Thing 

Since the establishment of the canonical JAK-STAT pathway, 
it has subsequently been replicated multiple times to establish 
the diversity seen within extant higher vertebrates. However, 
the pathway itself has not been replicated per se, but rather the 
individual components have undergone replication, creating an 
opportunity for the generation of new pathway combinations. 

Duplication events, either local or global, have largely facili- 
tated the expansion of gene families in vertebrates. At the global 
level, three whole genome duplications (WGDs) have molded 
vertebrate evolution: the so-called IR and 2R, which occurred 
before the divergence of lobe-finned and ray-finned fishes,"* and 
3R, which was limited to teleost fish (Fig. 3).'' While rapid loss 
of many duplicated genes occurred following WGD, some have 
been retained because of either division of function between the 
duplicated genes or the generation of new functions for one of the 
duplicates.'^'"^"'"^' We and others have shown that WGDs supple- 
mented by local duplications have generated the diversity of JAK- 
STAT pathway components seen in higher vertebrates. 

Following the establishment of the archetypal JAK (typified 
by fruit fly Hopscotch) prior to the divergence of protostomes 
and deuterostomes, there has been no increase in the number of 
JAK homologs until the emergence of vertebrates, with existing 
urochordates such as sea squirt {Ciona intestinalis) retaining a 
single JAK.*"^ The four JAKs in present-day tetrapods (e.g.. Homo 
sapiens) — ^JAKl, JAK2, JAK3 and TYK2 — were generated from 
IR and 2R via JAK1/TYK2 and JAK2/JAK3 intermediates, 
while 3R yielded the additional teleost-specific JAK2 paralogs, 
named jak2a and jak2b,'' such as in zebrafish [Danio rerio). 

The archetypal STAT (typified by fruit fly Marelle) was simi- 
larly finalized by the divergence of protostomes and deutero- 
stomes.''^ However, local duplication appears to have generated 
distinct duplicates in chordates prior to divergence of the urochor- 
dates, as observed in present-day sea-squirt (stat.a and stat.b).*"" 
We have argued that the proximity of many existing vertebrate 
57/1 r genes suggests these duplicates originally lay adjacent as a 
consequence of local duplication of the original STAT gene, with 
this STAT cluster being subsequently duplicated en bloc via IR 
and 2R. Subsequent gene loss has resulted in three discrete cop- 
ies of this cluster, STAT3-STAT5, STAT2-STAT6 and STATl- 
STAT4. Finally, 3R resulted in retained duplicates of both STATl 
(statl.a and statl.b) and STAT5 (stat5.1 and stat5.2), while local 
duplication has occurred along the mammalian lineage to gener- 
ate adjacent STAT5 genes (STAT5A and STATSB)." 



In contrast, the evolution of the SHP family appears to have 
been driven solely by WGD processes in concert with poor 
duplicate retention. The archetypal SHP (typified by fruit fly 
Corkscrew), established prior to the protosome and deuterosome 
split, was duplicated via IR and 2R, with subsequent gene loss 
yielding three members (SHPl, SHP2 and SHP3) in non-mam- 
malian tetrapods. The 3R WGD failed to generate any further 
expansion of SHP family members in teleosts, while SHP3 has 
been lost in the mammalian lineage.'' 

The common ancestor of protostomes and deuterostomes 
possessed three SOCS members: SOCS1/2/3/CISH, SOCS4/5 
and SOCS6/7.* A further duplication of SOCS6/7, to generate 
SOCS6 and SOCS7, occurred prior to the divergence of uro- 
chordates, as exemplified in present-day sea squirt. WGDs have 
mediated the on-going expansion of the family within verte- 
brates. SOCSl, SOCS2, SOCS3 and CISH were generated from 
the SOCS1/2/3/CISH precursor during IR and 2R, via SOCSl/ 
SOCS3 and SOCS2/CISH intermediates. Concurrent gene loss 
limited the expansion of the other genes to yield SOCS4 and 
SOCS5 from the SOCS4/5 precursor and no further expansion 
of SOCS6 and SOCS7. In contrast, the 3R event produced para- 
logs for SOCS3 (socs3a and socs3b), SOCS4 (socs4a and socs4b), 
SOCS5 (socs5a and socs5b) and CISH (cish.a and cish.b), as seen 
in the present-day teleost, D. rerio} 

Loss of duplicated genes is attributed to the generally redun- 
dant functions of the duplicate.'*''" Therefore, appropriate posi- 
tive selective pressure is required to maintain duplicates. We have 
hypothesized that the key selective force in the case of JAK-STAT 
pathway component retention has been the expansion of the 
upstream cytokine receptors. This occurred co-incidently with 
an increase in organism complexity, including the emergence 
of adaptive immunity, providing an opportunity for the dupli- 
cates to develop crucial, non-redundant functions that facilitated 
retention.'^ 

Conclusion 

There have been four critical phases in the evolution of the com- 
plex canonical JAK-STAT pathway: (1) the generation of the 
individual domains with discrete functions, such as the writ- 
ing, reading, and erasing of a phosphotyrosine code, or other 
accessory roles, most of which occurred early in eukaryotic evo- 
lution, (2) the arrangement of these domains into functional 
proteins, which happened gradually across many years, (3) the 
interaction of these proteins to form an active pathway, which 
occurred in bilateria, and (4) the subsequent proliferation of 
this pathway, which was largely driven by WGD during verte- 
brate evolution. This has collectively facilitated the generation 
of a complex JAK-STAT pathway able to orchestrate responses 
to appropriate extracellular signals to produce the coherent 
melody of transcriptional outcomes that underpins hematopoi- 
esis, immunity, growth and other essential functions in higher 
vertebrates. 
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